The generalization of the Regge- Wheeler equation for self-gravitating matter fields 
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It is shown that the dynamical evolution of perturbations on a static spacetime is governed by 
a standard pulsation equation for the extrinsic curvature tensor. The centerpiece of the pulsation 
equation is a wave operator whose spatial part is manifestly self-adjoint. In contrast to metric 
formulations, the curvature-based approach to gravitational perturbation theory generalizes in a 
natural way to self-gravitating matter fields. For a certain relevant subspace of perturbations the 
pulsation operator is symmetric with respect to a positive inner product and therefore allows spec- 
tral theory to be applied. In particular, this is the case for odd-parity perturbations of spherically 
symmetric background configurations. As an example, the pulsation equations for self-gravitating, 
non-Abelian gauge fields are explicitly shown to be symmetric in the gravitational, the Yang Mills, 
and the off-diagonal sector. (PACS numbers: 04.25.Nx, 04.40.-b, 04.70.Bw) 



In general relativity, several interesting and even some 
of the outstanding problems lie, at least partially, within 
the scope of gravitational perturbation theory. A promi- 
nent example in vacuum gravity is the collision of black 
holes, for which the close limit approximation to head-on 
collisions has proven to be in remarkable agreement with 
fully numerical simulations Q|, 0. In the presence of 
matter, the stability of neutron stars with respect to non- 
spherical perturbations ||^ , the study of critical phenom- 
ena in the weakly non-spherical gravitational collapse [Q , 
as well as stability and local uniqueness questions for 
hairy black-holes and self-gravitating solitons , g pro- 
vide further topic examples wherein perturbation theory 
finds a natural application. 

For vacuum gravity, perturbation theory is highly de- 
veloped and even for self-gravitating systems, important 
general properties of the perturbation equations were 
established. In particular, it was recently shown that 
a manifestly hyperbolic, gauge-invariant formulation of 
gravitational perturbation theory does exist for static 
background configurations 0, ^. However, only for spe- 
cial systems with linear matter models, such as vacuum 
gravity and the Einstein-Maxwell system M, formula- 
tions of perturbation theory are known for which, at least 
on a static background, the perturbations are described 
by a system of pulsation equations, that is, by a wave 
operator whose spatial part is (formally) self-adjoint. In 
view of the powerful results for self-adjoint elliptic oper- 
ators, such symmetric formulations of gravitational per- 
turbation theory would be most valuable for analytical as 
well as numerical investigations of the problems above- 
mentioned, in particular for spherically symmetric back- 
ground configurations. 

In this letter we show how to cast the equations gov- 
erning the dynamical perturbations of static background 
configurations into a system of pulsation equations. We 
do so by first following ideas introduced by Choquet- 
Bruhat et al. for full, non-linear general relativity 0. 
In a second step we show that the resulting linearized 
wave equation for the extrinsic curvature can be sym- 



metrized. The relevant inner product is positive definite 
on the subspace comprising all traceless perturbations. 
In contrast to the traditional metric approach to gravi- 
tational perturbation theory, the method introduced in 
this letter admits a natural generalization to non- vacuum 
perturbations of non-vacuum configurations. In order to 
illustrate why, in general, the metric approach fails to 
yield a self-adjoint system of wave equations, we start 
with a brief account of the gauge-invariant formalism in- 
troduced by Gerlach and Sengupta |1J. For simplicity 
we restrict ourselves to odd parity perturbations. 
The gauge-invariant metric approach - The odd par- 
ity perturbations of a spherically symmetric spacetime 
{M, g) with metric g = g-\-B?g are parametrized in terms 
of a scalar field k{x^) and a one-form h = ha{x'')dx'^ 
on the two-dimensional pseudo-Riemannian orbit space 
M = Af/S0(3); 

Sgab^O, SgAb ^ hbSA, 6gAB ^'^i^^ [aSb), (1) 

where Sa = [*AY^™-)a denote the transverse spheri- 
cal vector harmonics and 2V(^S'b) = V aSb + ^ bSa- 
[Lower-case Latin indices (a = 0,1) refer to coordinates 
on {M ,g), while capital Latin indices (^4 — 2, 3) refer to 
the standard coordinates on the metric sphere (5^,9)]. 
Under infinitesimal coordinate transformations induced 
by a vector field X^ , the metric perturbations transform 
according Xo k —> K-\-ip and h ^ h-\- R^d (i?~^(^) , where 
If is the scalar field parametrizing the odd parity trans- 
formations, X'^ — 0, X"^ = ipg^^SB. In terms of the 
manifestly gauge- invariant one- form H , 



H^h-R^di — ), 



(2) 



the gauge-invariant components of the Einstein tensor 
0, SG'Xb = -d^'HVt^ASB) for £ > I, and 



become (5GX 



SG'^^dx^ 



Sa 
2i?2 



dt 



RM 



H 
R^ 



XH 



(3) 



Here, X = {£ — 1){£ -\- 2) and d^ = ?d? denotes the co- 
differential operator with respect to g. 



The vacuum equations, i5G™" = 0, imply a wave equa- 
tion for a single scalar field: The Regge- Wheeler (RW) 
equation Ijl^l for the radial component of H is derived by 
using the integrability condition d* H = Q (li (. ^ 1) io 
eliminate the temporal component of H. Alternatively, 
one obtains a covariant derivation (which is not limited 
to stationary background configurations) by applying the 
operator ?d to Eq. (||). This yields 



□ + i.O|l)+^ 



* = 0, 



(4) 



where the scalar ^ = R^*d{H/R^) is gauge-invariant for 
all values of £, and Q denotes the d'Alembertian with 
respect to g. 

In the presence of matter fields the gauge-invariant 
Einstein tensor 6G™J' is still adapted for analyzing sta- 
tionary perturbations [pi. However, following the above 
lines, it turns out that the evolution equations for met- 
ric perturbations cannot be cast into a self-adjoint form. 
In fact, as the derivation of the Teukolsky equation sug- 
gests, it cannot be expected that the gravitational per- 
turbations are parametrized in terms of a single scalar 
amplitude, unless for vacuum (or electro-vacuum) per- 
turbations on a algebraically special background. It is, 
therefore, natural to look for a wave equation for the one- 
form H itself. However, the operator governing the evolu- 
tion of H is not hyperbolic, since the expression (||) com- 
prises only "one half" of the d'Alembertian. Hence, we 
are aiming at a description of gravitational perturbations 
for which the wave operator already appears "off-shell" . 
Laying stress on the staticity rather than the spherical 
symmetry of the background, the natural gauge- invariant 
perturbations turn out to be curvature quantities. 
The gauge-invariant curvature-based approach - Under 
infinitesimal coordinate transformations generated by a 
vector field (AT^) = (/, X'-), the metric perturbations of a 
static background, g = —a^dt^ -\-g, transform according 
to Sa ^ Sa + oi^iX^ + a/, Sgu — > Sgu -\- Xi — ct^ f,i and 
Sgij —* Sgij + Xi-j -\- Xj-i, where a dot denotes differenti- 
ation with respect to t, and covariant derivatives refer to 
the Riemannian background metric g. Hence, the tensor 
Sgij — 5gti-j — 6gtj:i is gauge- invariant under transforma- 
tions with / = 0. Considering odd parity perturbations 
on a spherically symmetric background from this point 
of view, the natural gauge- invariant amplitudes are given 
by the (lA) and the {AB) component of the above tensor. 
By virtue of Eqs. (|l|) these are proportional to 



hi - E?di 



i?2 



ht. 



(5) 



While the RW one-form H is suited to describe station- 
ary perturbations, the gauge-invariant variables (o) are 
adapted to their time-evolution. Note that the situa- 
tion is reminiscent of Maxwell theory, where the gauge- 
invariant quantity satisfying a wave equation is the elec- 



tric field, Ei = —Ai -\- <I>,i, which is constructed in a simi- 
lar way as the tensor Sgij — Sgu-j — Sgtj-i- Since the latter 
essentially is the perturbation of the extrinsic curvature 
Kij, our goal is to derive a symmetric wave equation 
for 6Kij. Within the curvature-based approach the RW 
equation has, as we shall show below, a two-component 
analogue which admits generalizations to self-gravitating 
field theories in a straightforward way. 
The symmetric wave equation for the extrinsic curvature 
- In the Arnowitt-Deser-Misner (ADM) formalism the 
metric reads 



-a^dt^ + gij{dx' + P'dt)idx^ + /3Mi), 



(6) 



where the shift /?* and the extrinsic curvature Kij = 
{2a)~'^ {[dt — Cp]g), ■ vanish on a static background. 
An efficient way to linearize the ADM equations is to 
use "vector-invariant" quantities, that is, perturbations 
which are invariant under the subset of transformations 
generated by vector fields (Xf^) = {f,X^) with / = 0. 
These quantities are 6Roi, SKij, and 



Sa — Cspo. 



i^spg). 



SR^ 



iCspR] 



ij ' 



(7) 



where the index "0" refers to the normal unit vector field 
do = a^^{dt — P^di). As the variation of the shift en- 
ters the ADM equations via the combinations (0) only, 
one may set (5/3 = for computations, provided that in 
the resulting expressions all perturbations are identified 
with their vector- invariant counterparts. The only gauge 
freedom is then parametrized by the function /. In par- 
ticular, one has 

[a^SK - Sa] -> [a^6K - 6a] - af + a^ {af-') .^ , (8) 

where now 5a denotes the vector-invariant lapse de- 
fined in (M). By virtue of this transformation law the 
scalar gauge freedom can be used to impose the har- 
monic gauge, a^5K — Sdt. Adopting this gauge and using 
vector-invariant quantities removes all gauge degrees of 
freedom, up to residual gauge transformations generated 
by {X^') = (/,0) with / = a{af ■')■,. The linearized 
ADM evolution equations then become 



SK, 



y - —^9^J, SR^j = 5R,j + — ^ - 5 — -^, (9) 



2a "■' " 'a a 

while the linearized constraints are given by 

<5Goo = \6R, 5Rq, = V^K,, - V,SK. (10) 

In order to obtain a wave equation for SKij we first 
follow the method introduced by Choquet-Bruhat et al. 
H, and then show how to cast their result into a sym- 
metric form. More precisely, one may proceed as follows: 
First, one differentiates the second evolution equation in 
Eqs. (PI) and uses the first one to obtain an expression 



for 6Rij in terms of 5Kij and covariant derivatives of 
6Kij up to second order. Second, one eliminates the 
variation of the lapse by adopting the harmonic gauge, 
6d = a^SK. Next, one gets rid of those second covariant 
derivatives of SKij which are not of the desired form by 
subtracting a[{SRoi);j + {SRoj)-i] from SRij and apply- 
ing the commutation law for covariant derivatives. The 
resulting expression, which still contains first covariant 
derivatives of SKij, is finally symmetrized by subtract- 
ing S[a-jdRQi + a-iSRoj]. The result is that the tensor 



A,: 



= a^dtSR, 



(a^SRoA,.- (a^SRt 



Ojj 



(11) 



yields a formally self-adjoint, hyperbolic operator acting 
on SKij. [The difference between the Ay defined in Eq. 



( |l l| ) and the flij introduced in Ref. 
metrizing terms.] One finds 



lies in the sym- 



Ay = aSK,^ + 2a^ R\^5Kj)i - 2a'' Ruk.SK'^' 

+2V(i [aa^] a5Kj)k - 2V, [aa^] af^SKu, (12) 



where \3SKi, 



\df — aW'^aWkCt]SKij. The complete 



operator in Eq. (|l2| ) acting on the linearized extrinsic 
curvature is symmetric with respect to the inner product 
induced by the Dc Witt metric 



Qtjkl ^-ik-jl _-tJ-kl 



(13) 



Hence, the inner product is positive definite only for per- 
turbations with vanishing trace, SK — 0. In particular, 
this is the case for odd parity perturbations of spherically 
symmetric background configurations. 

For a spherically symmetric background with metric 
g = a^dp^ + R^dfl^, the multipole expansions of the lin- 
earized vacuum equations, ApA = 0, Aab = 0, yield the 
pulsation equation [□ + P]u = and the (momentum) 
constraint (R'^a^^ui)' = ^/XRu2, where 



P = 



1 



i?2c 



R'^/a{a/R^)" + \a^ 2V\R^ia/Ry 
2%/Ai?2(a/i?)' RR" + Xa^ 



The gauge-invariant amplitudes u =^ (?ii,U2) are defined 
by SKpB = uiSb and SKab = ^ck^^A~^'2u2V(yiS's), 
the d'Alembertian refers to the background metric g = 
a2(— dt^ +dp2), and a prime denotes differentiation with 
respect to p. 

A convenient way to solve the initial value prob- 
lem is to impose initial conditions u{0,p) and u{0, p) 
subject to the constraint equation and its time deriva- 
tive. Having solved the wave equation for u{t,p), the 
metric perturbations are then obtained from the SRij 
equation by using the relations —21*2 = vXR^^Ht and 
— 2ui — aR?id{R~^ H) between the old and the new 
gauge-invariant amplitudes. 

Pulsation equations for self- gravitating matter fields - 
With respect to the ADM basis of one-forms the gauge 



potential A is parametrized by a scalar field $ and a 
one- form A (both Lie algebra valued) , 



A^^aAt + Ai {Ax' + (5' At) . 



(14) 



In order to linearize the YM equations we proceed in 
a similar way as in the gravitational case. While the 
fundamental gravitational object is the extrinsic curva- 
ture, Kij = ^Cg^Qij, the corresponding quantity in the 
YM case is the electric one-form E — —iggF. Here, 
F — AA + AaA'is the YM field strength and id„ denotes 
the inner derivative with respect to the normal unit vec- 
tor field do- For a static, purely magnetic background 
one has A — A, while $ and E vanish. The linearized 
electric field thus becomes 



aSE = D(a(5$) -SA + Cs(,A, 



(15) 



where D = d -t- [A, ] is the three-dimensional gauge co- 
variant exterior derivative. Since SA —f 6A + CxA under 
coordinate transformations, SE, (5$ and SA — CspA are 
vector-invariant quantities. Hence, the variation of the 
shift can be taken into account by using vector-invariant 
metric and matter perturbations. It then remains to 
linearize the YM equations with respect to the time- 
dependent metric g = —a^At^ + gijAx^'Ax^ and the YM 
potential A — ^aAt + A. As $ and the temporal deriva- 
tives of all background quantities vanish, the evolution 
equations become 



SA^I){a6^)-aSE, 5E ^ SB' [aF) , 
while the Gauss constraint is given by 



D^SE = 0, 



(16) 



(17) 



where, for Lie algebra valued p-forms, D = (— 1)p*D5. 
These equations hold for arbitrary variations of the shift, 
provided that all perturbations are identified with their 
vector-invariant counterparts. Moreover, SE is invariant 
under arbitrary gauge transformations generated by Lie 
algebra valued scalar fields Xi since 5^ ^ S^ + x and 
5A^SA + Bx- 

The linearized gravitational equations (||) , (^0|) and the 
YM equations (|l^), (|l7| ) have the same structure, where 
SKij and SE, as well as Scji^j and SA correspond to each 
other. In both cases one has to differentiate the sec- 
ond evolution equation and to use the first one to obtain 
hyperbolic equations for SKij and SE. Also in both situ- 
ations one needs to complement the resulting expressions 
with spatial derivatives of the constraints to obtain the 
desired wave operator. Performing this procedure for the 
YM equations (|6|), (^, and using the background equa- 
tion D (aF) = yields the linearized YM equation 



A^adt 



^metricD^ ("^) 



(18) 



The RHS describes the couphng to the metric perturba- 
tions Sa and 6gij, while the one- form A comprises the 
matter perturbations; 



A= a5E + a^*[SE,*F] 



(19) 



Both, 



where USE = [adf + Da^D^ + aD^aBalSE. 

the wave operator and the potential in Eq. (noh are 

symmetric with respect to the invariant inner product 

-AGjg'HiaceUE^^'>SE^'^'>\. 

It is remarkable that the linearized gravitational equa- 
tion 



Ay = 8ttG 



a^dtST,, 



^{aHT^,), 



{aHTo,) 



(20) 



and the linearized YM equation (|18|), with Ay and A 
given in Eqs. (|lj) and (|l9|), respectively, perfectly fit to- 
gether: All matter perturbations arising on the RHS of 
the gravitational equation (20) can be expressed in terms 
of SE by using the first YM evolution equation in ( |lq ) 
and its first spatial derivative. In a similar way, the RHS 
of the YM equation (nSf) involves only metric perturba- 
tions which can be expressed in terms of 6Kij by taking 
advantage of the first ADM evolution equation in (||). 
Hence, in the wave equation for SKij the matter pertur- 
bations arise via SE only, and vice- versa. Moreover, using 
the momentum constraint, the off-diagonal contributions 
arising from the coupling between the gravitational and 
the YM sector enter Eqs. ( p^ and ( pO| ) in a perfectly 
symmetric way. We thus end up with a single, formally 
self-adjoint wave equation describing both gravitational 
and matter field perturbations. 

The derivation of the RHS of Eqs. (|l|) and ^) will 
be given elsewhere |12| ]. Here we only present the result: 
The gravitational pulsation equation ( po|) becomes 



Ay = SGtraceja^ 



F^F\ 



^FkiF'^'SK, 



-n ■F'^ F'-' 

7kl -. 



SK, 



kl 



-g.ja^F^'Dk [aSEi] 



-Dk[a5Ei^^a'F^ 



i) 



with Kij according to Eq. (|1 
tion (111) is given by 



aSEkD. 



(* 



2 f=.fc 

a r 



3) 



} (21) 



The linearized YM equa- 



A, = ~8Ga^F\ trace [5EjF{\ + 2aD^ [a'^F'^^] SK^k 
-2a^F)V^ [aSK.k] + 2a^F,ja^5K (22) 

with Ai according to Eq. (19[). It is not hard to verify 
that the above equations form a symmetric system with 
respect to the inner products given above. The pertur- 
bations SKij and 5Ei have to satisfy the momentum con- 
straint and the YM Gauss constraint (1^, respectively, 
while 5gij and 5Ai are subject to the Hamilton constraint 
and Eqs. (|) and (|l|). 



For odd parity perturbations on a spherically sym- 
metric background, the multipole decomposition of Eqs. 
(p^), (p^, ( pi] ) and ( p^ ) yields a similar pulsation equa- 
tion as in the vacuum case, where u now comprises two 
gravitational and five YM amplitudes [for SU(2)]. The 
explicit form of this pulsation equation as well as a gauge- 
invariant formulation of all constraint equations will be 
given in Ref. Q. 

Conclusion - The dynamical perturbations of a static 
spacetime are governed by a system of standard pulsa- 
tion equations for the linearized extrinsic curvature. In 
the presence of matter fields, the curvature-based gravi- 
tational perturbations are naturally embedded in a larger 
system of standard pulsation equations, comprising both 
gravitational and matter perturbations. In particular, 
the dynamical perturbations of static EYM configura- 
tions are described by a wave operator whose spatial part 
is manifestly self-adjoint. As the relevant inner product 
in the gravitational sector is positive definite only for 
traceless perturbations, further investigations are needed 
to discuss, for instance, the dynamics of even parity per- 
turbations in a satisfactory way. For perturbations with 
vanishing trace, in particular for odd parity perturba- 
tions of spherically symmetric configurations, the form 
of the pulsation equation is optimal and enables one to 
apply the tools from spectral analysis. 
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